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Abstract: In this paper we present an analysis at NLO of the contribution from squark-squark 
production to the experimental signature 2j+l + l~ +]£t(+X) with opposite-sign same flavor leptons, 
taking into account decays and experimental cuts. We consider the case in which one squark decays 
directly into the lightest neutralino Xi an d the other one into the second lightest neutralino and 
subsequently into l + l~Xi via an intermediate slepton. On one hand we study effects of the NLO 
corrections on invariant mass distributions which can be used for future parameter determination. 
On the other hand we analyze the impact on predictions for cut-and-count searches using this 
experimental signature. 
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1. Introduction 

Supersymmetry (SUSY) remains one of the theoretically most appealing models for physics beyond 
the standard model (SM), and the light Higgs boson predicted by the Minimal Supersymmetric 
Standard Model (MSSM) is consistent with the recent observation of a Higgs-like bosonic resonance 
[1,2]. However, for a natural stabilization of the EW sector, direct signs of SUSY should hopefully 
be observable at the TeV scale. Until now there has not been any direct evidence for deviations from 
the SM and the experiments at the LHC have set new limits in various regions of the parameter 
space of the MSSM, investigating numerous final states and signatures [3-7]. 

Once a clear signal for physics beyond the Standard Model (BSM) is established, the character 
and parameters of the underlying model have to be determined. Only a precise and non-ambiguous 
determination of TeV scale parameters allows reconstruction of the underlying theory [8, 9] and, 
e.g. in case of SUSY, investigation of different breaking scenarios. Moreover, different BSM models 
result in similar signatures at the LHC, so even the determination of the general model would be a 
major challenge [10]. In recent years various techniques for these challenges have been developed. 
Many of these techniques rely on the occurrence of cascade decay chains where, due to an addi- 
tional symmetry, the lightest new particle (LSP in the SUSY case) is stable and, motivated by the 
possibility to provide a dark matter candidate, unobservable in the detectors. For a review of mass 
determination techniques see, e.g., [11]. 

In this paper we investigate NLO (S)QCD corrections to the "qll-chain", also known as "golden 
decay chain" , 



where a left-handed squark decays into a quark and a second lightest neutralino Xs>j which subse- 
quently decays via an intermediate slepton into a pair of opposite-sign same-flavour (OS-SF) leptons 
and a lightest neutralino Xi- For a systematic treatment we combine this decay chain at NLO with 



QL -> q X°2 -> q I* ILr "> 1 l± F X°i 



(1.1) 
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the production of a pair of squarks q L q' R , where the second squark decays directly into a quark and 



a xi, 



QR -> q Xi ■ (1-2) 



Thus we provide, for the first time, a fully differential description at NLO QCD of the contribution 
to the signature 2j + l + 1~ (OS-SF) + $t{+X) from the process 

PP^QLQR^Qxiq'^xi- (1-3) 

For heavy squarks and gluinos, as suggested by LHC exclusion limits, the considered squark-squark 
production is the dominant process among all coloured SUSY production channels, see for exam- 
ple [12]. In our analysis we assume xi an d X2 t° be mainly bino- and wino-like, as they appear in 
large parameter regions of models with unified gaugino masses at the GUT scale. In such scenarios 
the other decays, qx — > XiQ an d Qr X2Q1 are highly suppressed for squarks of the first and second 
generation. Furthermore, in the benchmark scenarios we consider, all relevant squarks are lighter 
than the gluino g, so squarks can exclusively decay into neutralinos and charginos. 



The decay chain, eq. (1.1), was introduced in [13,14] and studied in many subsequent works [15- 
32]. These analyses showed that measurements of resulting invariant mass distribution endpoints 
and shapes can be exploited to determine the masses of the intermediate SUSY particles. Shapes of 
various relevant invariant mass distributions have first been calculated analytically at LO in [18] and 
they might be important to resolve ambiguities in mass measurements from kinematic endpoints 
[33]. In figure ^, for illustrative purposes, normalized invariant mass distributions in rriju, wijz(high) 
and TOj7(io W ) are displayed for two benchmark points (with details given in the corresponding sections 
of this paper). At this point we only want to emphasize that such shapes are very sensitive to the 
model parameters; on the other hand, they can be distorted by higher order corrections, affecting 
thus the accuracy of endpoint determination techniques. 




[GeV] [GeV] 

Figure 1: Shapes of the rriju, m j; ( high ) and mji(i ovr ) distributions for SPSla and 10.1.6 at LO. 



Such invariant mass distributions, due to correlations between particles in the decay chain, can 
be used also for spin measurements of the intermediate sparticles [34]. In this way, e.g., a SUSY 
model can be discriminated from an Universal Extra Dimension (UED) model, see e.g., [35, 36] 
and [37-47]. Additionally, besides observables based on invariant mass distributions, the signature 
2j + / + ?~(OS-SF) + ]£t{+X) can be used for SUSY searches at the LHC. Basic kinematical cuts 
can reduce SM backgrounds significantly and already now stringent bounds on relevant parameter 
regions have been obtained [3,6]. Furthermore, resulting rates can also be used for parameter de- 
termination within a global fit [48] . 
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In all these analyses a detailed understanding of theoretical uncertainties and effects from 
higher-order contributions is necessary, requiring precise predictions at a fully differential level 
including NLO corrections in production and decay. In [20] the "qll-chain" was investigated at 
NLO QCD in the squark rest frame. There, real gluon radiation contributions are given in a fully 
analytical form and leading soft- and collinear gluon contributions are resummed. Furthermore, 
an LHC analysis using LO production matrix elements is presented. In their numerical analysis 
the authors of [20] concentrate on the effect on distributions sensitive to spin correlations. In the 
present paper we continue the analysis of higher-order corrections of the "qll-chain" by combining 
the production of squark pairs and their subsequent decays at the NLO level in a consistent way. 
This allows a systematic study of the important jet combinatorics issues in the presence of more 
than one final-state jet, emerging from both real radiation at NLO and the combination of decay and 
production [29]. Futhermorc, we investigate, besides the effects on spin determination, the impact 
of such corrections on distributions important for mass determination and on inclusive event rates 
after experimental cuts. 

Many details of the calculation, not explicitly discussed here, can be found in [49] where an 
analogue analysis at NLO QCD of production of squark pairs directly decaying into two lightest 
neutralinos has been performed. Also another similar calculation for stop production and direct 
decay was recently performed in [50]. Conversely, calculations of higher-order contributions to 
either production [12, 51-72] or decay [73-85] of squarks (and gluinos) are manifold and include 
electroweak NLO corrections as well as QCD corrections beyond NLO. The study of the signature 
2j + l + l~ (OS-SF) + ]£t(+X) presented here includes the contribution from squark-squark produc- 
tion, which constitutes the dominant production channel in the case of heavy squarks and gluinos, 
as already pointed out. This work should be understood as a first step towards a calculation in- 
cluding all channels and also off-shell and non-factorizable corrections. 

The outline of this paper is the following. In section |2| we briefly explain the method, based 
on the previous study [49], to combine production and decay at NLO. In section pi we describe the 
calculation of the NLO corrections to just the decay chain, eq. ([hi]). In section^ we present the 
numerical results. This section |] is dived in three subsections: the discussion of the considered 
benchmark scenarios, the analysis of NLO corrections to the isolated decay chain, and finally the 
discussion, for combined production and decay, of NLO effects on invariant mass distributions and 
also on event rates after kinematical cuts. Conclusions and a possible outlook are given in ||. 



2. Method 

We investigate the production of (Ilq'r (QlQr) pahs induced by proton-proton collisions, where the 
right-handed squark directly decays into the lightest neutralino Xi an d the left-handed one into 
a x® an d subsequently into l + l~x\ via an intermediate slepton l^/ R - This process results in the 
signature 2j + l + l~ + J£t{+X) and is illustrated in figure |^. 

At LO the only partonic subprocess that contributes to a given intermediate <7l<?/j (QlQr) 
configuration arises from a quark (anti-quark) pair qq' (qq 1 ) in the initial state. In our calculation 
we include contributions from all (s)quark flavors of the first two generations. We will perform 
the following discussion without referring to the charge-conjugate subprocesses; in the final results, 
however, we include the charge-conjugate processes explicitly. Indeed, due to chirality-dependent 
interactions, the decay of an anti-squark has to be treated independently from the corresponding 
squark decay. 

Using the same theoretical framework as applied in [49] for squark-squark production and direct 
decays into lightest neutralinos, we include the NLO factorizable corrections in the narrow-width- 
approximation (NWA) to the process defined in eq. (|l.3| ). They correspond, as illustrated in figure 
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Figure 2: General structure of factorizable NLO QCD corrections to the given process. 



to the separate sets of corrections for the production and for the decays with squarks treated on- 
shell. For any flavor configuration, a systematic expansion of the differential cross section in the 
strong coupling a s yields 

qxWi + rx a i(+x)) = 



da< NWA^ pp ^ 
1 



(2.1) 



r (0) r (o) 

fc q' R 



da 
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In eq. (|2.l|), r^ ' and T-?' denote the LO total widths of the two squarks. The terms dT^ ~ , , , 

' c — 9i 9 H M qL^qx"l + l~ 



and (irl^ , ~ are the LO differential distributions for the decay chain and for the direct decay 
q R ~^q X\ 

into the lightest neutralino, respectively, boosted to the moving frames of qi, and q' R . The other 



ingredient, da. 



(0) 

pp^rqLq' R 



is the LO hadronic differential production cross section. The corresponding 



NLO corrections to the aforementioned quantities are indicated by an apex (1). The second line 
of eq. (^j]) contains a global factor with NLO contributions to the decays widths; the third line 
includes, as the first term, the corrections to the decay chain and, as second term, the corrections 
to the decay into the lightest neutralino; the last line represents the corrections to the squark- 



squark production cross section. The terms in eq. (2.1) refering to the LO and NLO contributions 
to squark-squark production and to the direct squark decays into the lightest neutralino were 
calculated and discussed in [49]; the treatment of the decay chain is addressed in section 0. On 



the basis of eq. (2.1) intermediate events for production and decay are combined in analogy to the 
strategy explained in [49]. 



3. Calculation of the squark decay chain 

The structure of the decay chain is illustrated in figure ||. With lf L and I J we indicate, respectively 
the lepton emerging from the \2 decay, the near lepton, and the lepton emerging from the slepton 
decay, the far lepton. Experimentally these leptons are indistinguishable. 

Figure ^ represents, in a compact notation, the four Feynman diagrams contributing to the tree 
level amplitude for the decay gx — > ql + l~\\. They correspond to the two cases qx — > <ll£lJXi an d 
qL Q^n^Xi with a left- or right-handed intermediate (anti) slepton. 
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For non-degenerate left- and right-handed 
sleptons, as in the scenarios investigated in this 
paper, in NWA the structure of the squared am- 
plitude of this decay chain becomes much sim- 
pler. In the limit ^ — > for the sleptons and for 
X2 , the interferences between different diagrams 
vanish and the phase-space of the intermediate 
particles can be treated on-shell. In this way LO 
and NLO contributions to the differential distri- 
bution for the decay chain in Figure pi can be written as follows 




Figure 3: Structure of the decay chain. NLO QCD 



corrections involve only the first step qL 



QX2- 



(04) 



dr (0,D 



dY 

I 

er = ±l/2 



X2 



+ dr~ 



(3.1) 



The index a represents the helicity of x°>- We keep, at LO and NLO, the full helicity dependence. 
Technically, we use, both at LO and at NLO, the matrix elements for the entire chain and set 
consistently the different on-shell conditions according to the intermediate states for the various 
contributions. In this way, the sum over the helicity states is automatically performed and off-shell 
effects can be switched on easily in a further study. The two terms contained in the square brackets 



of eq. (3.1) correspond to the two different charge configurations for near and far leptons in the 
decay of Both configurations get contributions from the left- and the right-handed slepton. 

As mentioned above, we concentrate on scenarios where the X2 * s dominantly wino-like. Thus, 
in the considered decay chain the coupling to a right-handed slepton, Ir, is heavily suppressed 
compared to the corresponding decay chain via a II- Only if the decay into a left-handed slepton 
is kinematically forbidden, mj > m^o > mj , the decay via a Ir can contribute substantially. 
Basically also the decay of the X2 a OS-SF lepton pair via a Z boson contributes when kine- 
matically allowed; these effects are, however, numerically not significant for the benchmark points 
considered in this paper. 

is performed following the procedure explained in [49] . drl^__^,_ 



The calculation of dT. 



includes contributions from loop corrections and from real gluon radiation. Consequently, in 
collinear and infrared safe regions, our NLO corrections are fully differential in the momentum 
of the gluon emitted. The rest of the electroweak decay chain is not affected by NLO QCD correc- 
tions. 



4. Phenomenological results 

In this section we present the numerical results of our calculation. First, we specify input parameters 



and relevant observables. Second, NLO corrections to the decay chain, eq. (1.1), not combined with 
the production, are investigated in the squark rest frame. Third, we present results for the decay 
chain combined, according to section^, with o_l^r production. 

4.1 Parameters and observables 

Standard Model input parameters are chosen according to [86]. As PDF set we use CTEQ6.6 [87] 
with the associated a^ ls (fiR) at NLO. Renormalization scale /ir and factorization scale (if are both 
set to the average mass of all light-flavor squarks, fj, = (mf = fJ>R = TUq. 

Numerical results in this paper are presented for two representative benchmark scenarios and 
the LHC with a center of mass energy of y/S = 14 TeV. We choose the two CMSSM scenarios SPSla 
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u L 


u R 


d L 


dR 


9 


h 


Ir 


x° 2 


X? 


SPSla 


563.6 


546.7 


569.0 


546.6 


608.5 


202.4 


144.1 


180.2 


97.0 


10.1.6 


1531.7 


1472.2 


1533.6 


1466.1 


1672.1 


536.6 


340.6 


592.4 


313.3 



Table 1: On-shell masses of the first generation squarks and sleptons, the gluino and the lightest and 
second lightest neutralino within the different SUSY scenarios considered. All masses are given in GeV. 



BR (%) 


QR X° 


QR -> X.2 


Ql — >■ Xi 


-t x.2 


X°2 it 


x°2 it 


X°2^Z 


SPSla 


98.5 


1.0 


1.5 


31.2 




13.1 




10.1.6 


99.8 


0.03 


1.5 


32.1 


28.4 


0.2 


0.2 



Table 2: Branching ratios for the decay of squarks into Xi an d X.2 an d for the decay of a X2 into right- 
and left-handed sleptons. Squarks and leptons of the first two families are considered, where branchings 
into second and first generation sleptons and their charge-conjugate contributions are summed. 



and 10.1.6 defined in [88,89]. The scenario SPSla has already been excluded by searches at LHC. 
However, it still serves as viable benchmark scenarios, where many detailed studies are available 
in the literature. The scenario 10.1.6 is still viable and can be tested in the near future. The low 
energy spectrum for both scenarios has been obtained with the program SOFTSUSY [90]. Sparticlc 
on-shell masses relevant for our analysis are listed in table |l|. Non-vanishing Yukawa corrections 
implemented in SOFTSUSY result in small mass splittings between first and second generation squarks 
and sleptons. We verified that, for the CMSSM scenarios SPSla and 10.1.6, the phenomenological 
effects originating from these small mass splittings are negligible in the study presented here. Thus, 
we set all second-generation masses equal to their first-generation counterparts. 

In both scenarios the gluino is heavier then all light flavor squarks. Thus, all these squarks decay 
exclusively into charginos and neutralinos. In table ^ corresponding branching ratios, calculated 
with SDECAY [91], are listed 1 . The right-handed squarks decay dominant ly directly into the bino-like 
X?)the left-handed squarks into the x° an d the lighter chargino, Xi ■ 

In table ||we also list branching ratios for the second lightest neutralino into light flavor sleptons. 
Branching ratios into first and second-generation sleptons are identical, and in table ^| we sum those 
contributions. For benchmark point SPSla only the right-handed Ir is lighter then the x.2- Thus, 
next to the decay into a r-slepton, this is the only available two-body decay. In our numerical 
analysis of SPSla, both for the decay chain alone and combined with the production, only the 
decay via a right-handed slepton is considered. In contrast, for 10.1.6 both sleptons are lighter than 
the j(2- Due to its wino-like nature the x.2 here decays dominantly into the left-handed II despite 
the smaller mass of the Ir. For simplification in our numerical analysis of 10.1.6 only the decay via 
a left-handed slepton is considered. As can be seen from table the contribution from the decay 
into a Ir (and also into Xi^) can be neglected safely. 

In all numerical results presented in the following we employ the anti-fey jet clustering algorithm 
with a jet radius of R = 0.4 implemented in Fast Jet 3.0.2 [92]. Furthermore we define a jet to 
fulfil the cut conditions 

p£>20GeV, |77 i( | < 2.8. (4.1) 

Thus, we arrive at an experimentally well defined result. When analyzing combined production and 

1 In tabic |i|wc list the average of the value of the branching ratios for up and down type squarks, that, however, 
differ at most by ~ 1%. Differences between branching ratios at LO and NLO for squark decays are negligible (less 
than per mill) for the considered scenarios and so not shown. 
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decay in section 4.3, the following realistic experimental cuts are applied: 

pi > 150 GeV , P J 2 > 100 GeV , \ m | < 2.5 , 

pf la > 20 GeV (OS-SF) , fi T > 100 GeV , (4.2) 

where we implicitly require the two leptons to have opposite charge and same flavor (OS-SF). 
Such cuts efficiently reduce SM backgrounds [3,6,16]. Furthermore, we assume that contributions 
from leptonic decays of r-leptons (from SM processes or the corresponding signal decay chain with 
an intermediate tau sleptons), charginos and bosons are removed in the standard way by 
subtracting events with opposite-sign different-flavor lepton pairs (OS-DF), see e.g. [11,13,15-17]. 
Thus these cuts help to isolate the decay chain under consideration. 

4.2 Squark decay chain 



Here we want to investigate NLO corrections to the isolated decay chain eq. (1.1) evaluated in 
the squark rest frame. In figure ^ we show various invariant mass distributions in the final state 
leptons and jet(s) for the benchmark scenario SPSla, LO distributions are shown in black and NLO 
distributions in red. As explained in section [j], shapes of such distributions are very important for 
the determination of masses and spins of sparticles. In order to highlight NLO corrections purely in 



the shapes, here and in section 4.3 we show all distributions, both at LO and at NLO, normalized 
to unity. 

Two kinds of combinatorial problems arise in studying invariant mass distributions involving 
the final state leptons and jet(s). First, as already mentioned in section [TJ, from an experimental 
point of view we cannot distinguish between the near and the far lepton on an event-by-event 
basis. This is a well known problem and many solutions have been suggested in the literature, 
see e.g. [11]. Second, it is not obvious which jet to choose to built the desired invariant mass 
distributions. Considering only the isolated decay chain starting with q~L, at LO only the jet from 
the squark decay is present in the final state. But at higher orders, due to real gluon radiation, 
further jets can be present. Here we always choose the hardest available jet to build the invariant 
mass distributions, as done in [20]. 

In the upper left/right part of figure ^ we show (unobservable) distributions in the invariant 
mass of the hardest jet and the negatively /positively charged near/far lepton from the decay chain 
ul — » jlnl~fXi- Here and in the remainder of this section we do not apply any cuts, but the 



jet definition cuts in eq. (4.1). In the center left/right part, on the other hand, we show (again 
unobservable) distributions in the invariant mass of the hardest jet and the negatively/positively 
charged far /near lepton from the decay chain with the opposite charges for far and near leptons, 
ul — > jl£lJXv Finally, in the lower part of figure |] we show, in some sense, the sums of the 
two previous contributions. These distributions are in principle experimentally observable (after 
combination with the corresponding production process). The lower left /right panel shows the 
invariant mass of the hardest jet and the negatively/positively charged lepton summed over near 
and far contributions (rriji- and rriji+). 

In the case of the decay of a left-handed anti-squark, all the distributions introduced so far 
are equal to the charge-conjugate ones of the corresponding squark, e.g, the m -,+ distribution from 
an u* L decay chain is equal to the m . ; - distribution from an decay chain. Hence, the analogue 
of figure [| for q* L would present the shapes of the distributions of the left column exchanged with 
the ones of the right column. These differences between squarks and anti-squarks obviously do not 
appear for quantities that are inclusive in the different charges of the leptons. 

In all of the plots in figure || NLO corrections tend to shift the distributions to smaller invariant 
masses, however, locations of endpoints are unaffected. Kinematical edges in the NLO predictions 
are rounded off compared to LO, still, overall shapes of the considered contributions seem to be un- 
altered. Results and distributions of the same type, as already stated, have been calculated in [20]. 
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Figure 4: Normalized differential distributions for SPSla in rriji,, rriji n for the two unobservable decays 
ul — > jlnlfXi (upper two) and ul — > jlnlJXi (central two) and in rriji+ and m^- (lower two) where 
contributions from the two decays are summed. LO predictions are shown in black, NLO in red. 



In their numerical evaluation the slightly different parameter point SPSla' was investigated result- 
ing in LO shapes somewhat different to those presented here. Qualitatively the NLO corrections 
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Figure 5: LO (black) and NLO (red) normalized differential distributions for SPSla in mj((i ow ), Jn^high), 
rriju and m j;; (thresh) (from top left to bottom right) for the decay chain cf^ — ¥ q ~ *■ 1 1 — *■ 9 ' ^ Xi- 

shown in [20] for SPSla' and ours for SPSla agree. Moreover, we also investigated distributions 
and corrections for SPSla' and found agreement with [20] 2 . 

In figure ||we look at a different set of invariant mass distributions again for SPSla: f7ij;n w)) 
TO jj(high)j Trijii and TOj-;/(thresh)- For TOjj(i ow ) and TO^high) w e select on an event-by-event basis the 
smaller and higher invariant mass between one of the leptons and the hardest jet. rriju is the 
invariant mass between the hardest jet and the two leptons and WjM(thresh) is the same distribution 
where an additional constraint on the invariant mass of the two leptons, — < to/;, is applied. 
Here, TO™ ax is the well measurable endpoint of the dilepton invariant mass distribution 3 . All these 
invariant mass distributions are in principle experimentally measurable and have extensively been 
discussed in the literature [11,13,16,17] (and references therein). From a measurement of their 
upper (and in case of TO-^^thresh) the lower) endpoints one might be able to extract relations for 
the masses of the intermediate sparticles. These relations often show ambiguities and particularly 
measurements like the threshold of TOjz;( t h re sh) might help to resolve these. Alternatively, shapes 

2 In reference [20] a different jet algorithm is used. Results for SPSla' presented there for y c = 0.002 agree best 
with our results obtained using the anti-fey jet clustering algorithm. 

3 In our numerical analysis we use the theoretical endpoints m™ ax = 80.0, 203.8 GeV for SPSla and 10.1.6. 
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Figure 6: LO (black) and NLO (red) normalized differential distributions for 10.1.6 in mji(i ow ), m.,;(high), 
rrijii and m.,;; (thresh) (from top left to bottom right) for the decay chain <jl — > q X2 — > 1 I ^ 1 I F Xi- 



of the presented invariant mass distributions might help to overcome these difficulties. Let us now 
look at the NLO corrections to these invariant mass distributions. Overall, again, distributions 
are shifted to smaller invariant masses. Also, upper kinks of rriju and rriju (thresh) are rounded 
off. These shifts might result in a slightly lowered accuracy in the measurement of the upper 



endpoints. Furthermore, the threshold of the m 



j/£(thrcsh) 



is diluted due to NLO corrections 4 . A 



precise measurement of this observable seems to be questionable. 

The same set of invariant mass distributions, TOjz(i ow ) 5 m ji(high), m jii and mjiitthresh), is shown 
in figure ^| for the parameter point 10.1.6 and the corresponding decay chain involving a II- For the 
main part, again, NLO corrections shift the differential distributions to smaller invariant masses 
and round off the upper kinks. Particularly for mjm w) this might result in a smaller possible 
experimental accuracy for determining the upper endpoint. Apart from the rounding off of the 
kinks, general shapes of the distributions are mostly unaltered by NLO corrections. Also, a possible 
dilution of the lower endpoint of mj^thresh) d ue to NLO corrections seems to be less severe for 10.1.6 
compared to SPSla. 



4 The theoretical lower endpoint of the m 



jll (thresh) 



JU(thresh) 



: 437.1 for 10.1.6 [15]. 



distribution is given by ^^j^hresh) 



215.4 for SPSla and 
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4.3 Combination of production and decay 

Now we want to investigate the combined process, eq. ([l^), where cuts defined in eq. (fh^) have been 



applied. In section 4.3.1 we first examine the impact on various differential distributions important 



for parameter determination. Afterwards, in section |4.3.2| we investigate the impact of NLO QCD 
corrections on inclusive OS-SF dilepton observables and thus on searches currently performed at 
the LHC. 

4.3.1 Invariant mass distributions 

When combining production and decay already at LO a combinatorial problem arises when looking 



at invariant mass distributions. As already mentioned in section 4.2 it is not clear which jet to 



choose for building the different invariant masses. Just choosing the hardest jet as was done in 



section 4.2 does not seem to be sensible, since the jet from the second decay, present already at LO, 
is often the hardest one. This is a well known problem in the application of the endpoint methods for 
mass determination and various methods have been developed to reduce this ambiguity. The easiest 
method is to always use the jet which gives e.g. the smaller m.ju value. In this way one improves 
the measurements of the upper endpoints without losing statistics, however, shapes are heavily 
distorted already at LO. In principle there are advanced techniques to solve this problem, amongst 
others, full kinematic event reconstruction [93] or hemisphere techniques [94-96]. However, these 
techniques are quite involved, parameter point dependent and not generically applicable. Here, 
we apply consistency cuts, also discussed in [17] to reduce the impact of the jet combinatorics 
ambiguity. Applying such consistency cuts means, we accept only events where one jet ji out of the 
two hardest jets ji,jk gives an invariant mass smaller than m^ z ax and the other jet jk an invariant 

K 

jit 

™>hU < ™fi? x < m jkU . (4.3) 



mass larger than m m& 



Now ji will be used in the following to build the invariant mass distributions. This technique is 
very efficient in reducing the jet combinatorics ambiguity, however, event rates are also reduced 



(see section 4.3.2 and particularly table 



In figure ^ we show for SPSla the same invariant mass distributions as already shown in figure 



Here, production and decays are combined at NLO, cuts of eq. (4.2) and further consistency cuts 
are applied. Due to the NLO corrections distributions are in general shifted to smaller invariant 
masses. Comparing just LO predictions in figure [j] and figure || particularly rriju and mj;;( t h rcs h) 
show a slightly different behavior introduced by the consistency cuts: the plateau is less prominent 
in figure |?]. Here, again we observe a dilution of the threshold in the Wj;;( t h r0 sh) distribution at 
NLO. Similar observations can be made looking at figure @ (and comparing with figure |) for the 
combined results of parameter point 10.1.6. Overall, changes to the shapes are moderate and, 
concerning the measurements of the upper endpoints only for TOj;(i ow ) an experimentally detectable 



effect is expected. Here wc want to note one thing: the consistency cuts, eq. (L3), are based only 
on rriju and this is why, already at LO, we observe contributions also beyond the theoretical upper 
endpoint. This effect is enhanced at NLO. 

From the discussion above usefulness of the threshold of »7ijzz(thresh) seems questionable. Ad- 
ditionally, a measurement of a lower endpoint is always subject to large experimental backgrounds 
[19,25]. As this threshold was introduced to solve ambiguities in the mass determination due to 
the near-far indistinguishability, new techniques for this purpose have been invented. In [25] the 



5 From an experimental point of view the endpoint m^f x is assumed to be measured in a first step where for example 
always the jet is chosen yielding the smaller rriju. Here, we use the theoretical endpoints m^f* = 450.6, 1147.7 GeV 
for SPSla and 10.1.6. 
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Figure 7: LO (black) and NLO (red) normalized differential distributions for SPSla in rriji(i ow ), w,j((hi g h), 
rriju and rriju (thresh) (from top left to bottom right) for combined production and decay. 



authors argue that all invariant mass distributions used for mass determination from the given 
decay chain should be built symmetrically under the interchange Z noar O / far . In this spirit they 
introduce a new set of invariant mass distributions, = m^i U > m2 i(d) = \ m jl ~ m2 jif\i 

TftjUa} = m 'ji n + m ji f an d m %(p) = m jln ' m jjj- Here we study the impact of the NLO QCD correc- 
tions on this class of distributions. In figure [|wc show the normalized LO and NLO distributions in 
m jl(u)' m jl(d) ,m jl(s) an< ^ m ji(p) a S ams t a quadratic scale. Shapes of these distributions are slightly 
changed due to NLO corrections, however, the possibility of measuring their upper endpoints (both 
endpoints in the case of rnju u -\) seems to be unaffected. 

Besides for mass determination, the given decay chain can also be used for spin determination 
or, more precisely, for spin distinction. As pointed out in [34] and many subsequent works, the 
asymmetry between the rriji+ and rriji- distributions defined as 

A _ d<j/dm jl+ - da/drriji- ^ 
da /drriji+ + da/drriji- 

can help to discriminate between a SUSY model and other models like Universal Extra Dimensions 
(UED) with a similar decay chain where the intermediate particles have different spins [35]. In 
figure we show LO (black) and NLO (red) predictions for this asymmetry for SPSla (left) and 
10.1.6. Again, cuts of eq. ( |4. 2[ ) and consistency cuts have been applied. The potential for spin 
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determination and/or model discrimination seems to be unaltered by NLO QCD corrections. For 
SPSla at NLO there is a contribution beyond the upper endpoint not present for this observable 
at LO. This NLO contribution passes the consistency cuts which are based on the endpoint of the 
rriju distribution. However, expected event rates in this region of the asymmetry distribution are 
experimentally negligible. 



7^ ' ' ' I 1 n | I I r-p i r-p r 




200 400 600 800 1000 



m ii ( iow) [GeV] 




200 400 600 800 1000 1200 1400 
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Figure 8: LO (black) and NLO (red) normalized differential distributions for 10.1.6 in m j; ( low ), TrijiQaM, 
rriju and Trijim^resh) (from top left to bottom right) for combined production and decay. 
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Figure 9: LO (black) and NLO (red) normalized differential distributions for 10.1.6 in mji( u ), rriji^), mji(s) 
and mji( p ) (from top left to bottom right) for combined production and decay shown with a quadratic scale. 
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Figure 10: LO (black) and NLO (red) normalized differential distributions for SPSla (left) and 10.1.6 
(right) in the asymmetry A for combined production and decay. 
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»r(0),cons. cuts 




cuts 

^23 + 21 + ^7- 






SPSla 


38.2 fb 


23.0fb 


1.36 


1.23 


1.34 


1.28 


10.1.6 


0.628 fb 


0.243 fb 


1.46 


1.39 


1.44 


1.41 



Table 3: LO iV (0) and NLO N {0+1) cross section predictions and K-factors Kn for the two benchmark 
scenarios SPSla, 10.1.6 at a center of mass energy vS = 14 TeV where the cuts of eq. (4.2) are applied. 
For comparison we also list the inclusive NLO production K-factor K pp ^qqi and K pp ^ f ^ L ^ . 



4.3.2 Inclusive observables 



As discusssed in the previous sections, the experimental signature of two jets, two OS-SF leptons 
and missing energy provides essential information for future SUSY parameter determination. Ad- 
ditonally this signature can also be used for searches for supersymmetry. At the LHC, using this 
signature, cut-and-count searches have been performed in the analyses of the 7 and 8 TeV runs [3,6] 
and will be performed also at 14 TeV. Consequently, precise calculations of inclusive cross sections 
for the specific signal region used in the cut-and-count searches are necessary. 

First, we discuss the differences between predictions in our approximation and rescaling LO 
with a flat K-factor from NLO corrections to squark-squark production without including decays 



and cuts. Second, we look at the impact of the additional consistency cuts defined in eq. (4.3). In 



table |] various integrated quantities at 14 TeV for the parameter points SPSla and 10.1.6 are listed. 
Starting from the first column on the left we display: the total cross section at LO N^ +21+ ^ t in 



the signal region defined by the cuts of eq. (4.2), the total cross section at LO after consistency cuts 
-^2i+2°+^T UtS ' together with the corresponding K-factors KN 2j+2l+f , T and ^™ 2 "+ 2 ^f ■ Furthermore, 
we list the K-factors for the production only including all Qlq'h channels Kpp^q L q' R an d including 
also all other chirality configurations K pp ^qqi . 



The difference between the K-factors including the cuts defining the signal region, K 



and the K-factors for production of qlQr pairs, Kpp^q Lq - R i is small, namely 2%. This difference 
increases to 8% for SPSla and 5% for 10.1.6 if the K-factor for just the production includes also 
the other chirality configurations, K pp ^qqi . Thus, for the scenarios analyzed here, NLO corrections 
can be safely approximated rescaling LO predictions with the K-factor obtained for the production 
part, provided that only the contributing chirality configurations are included in the calculation of 
the K-factor. This feature, however, cannot easily be generalized; for example, as can be seen from 
table |^, applying consistency cuts increases the differences between these two approximations. 

Consistency cuts are designed for the study of distributions relevant for parameter determina- 
tion, as the ones discussed in the previous section. However, they decrease the cross sections, as 
can be seen comparing N^ + 2i+m T an d N^^i+Kt^ •> without adding any obvious benefit in the 
context of searches. On the other hand, so far, we did not discuss the different normalization of the 
LO and NLO distributions shown in the previous sections. The values of ^jv 2 n + 2 ™p are exactly 
the ratios between the normalization of the LO and NLO results. These values are smaller than 
the K-factors obtained without consistency cuts. As discussed, such cuts are ideated to solve the 
jet combinatorical problem. Doing so, they also reduce (positive) contributions from real radiation 
of a gluon or a quark at NLO and consequently reduce the resulting K-factors. Thus, in contrast 
to the case without consistency cuts applied, fully differential factorizable corrections have to be 
taken into account for a precise estimation of NLO effects. 
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5. Conclusions 

In this paper we provided, for the first time, an analysis at NLO of the contribution from squark- 
squark production to the experimental signature 2j + + ]£t{+X), taking into account decays 
and experimental cuts. We focused on the impact of NLO corrections on invariant mass distributions 
that can be used, in case of discovery of supersymmetric particles, for parameter determination. 
We observe that general shapes, besides smoothing of edges and kinks and a shift towards smaller 
invariant masses, are not strongly altered. This seems to be an universal behaviour despite the 
strong dependence of shapes on the parameter region. 

We also analyzed the impact of NLO corrections including decays on the predictions for cut- 
and-count strategies used in discovery searches. The predictions depend on the cuts applied and 
in general can be different from the result obtained rescaling LO predictions with flat K-factors 
obtained from the cross-section of just production without decays and cuts included. However, in 
particular cases, results obtained in this approximation and using our calculation can be in very 
good agreement, provided that only the contributing chirality configurations are included in the 
calculation of the flat K-factor. 

In our framework, results can be easily extended including off-shell effects in the electroweak de- 
cay chain. Also different kinds of experimental signatures emerging from different decay chains can 
be analyzed easily. In general an analogue calculation of such signatures from the other production 
channels of squarks and gluinos is desirable. 
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